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1 . Introduction 


In Part 1 (SH4.2-4) we have quantitatively shown that our drift 
model can simultaneously fit the observed 1965-66 and 1977 electron and 
proton spectra, in contrast with spherically symmetric models. We 
interpret this result as evidence of a charge-sign dependent effect due 
to particle drift. Using the same set of parameters, we could also 
simulate a shift in the phase and amplitude of the diurnal variation in 
the cosmic-ray (CR) intensity by reversing the polarity of the 
interplanetary magnetic field (IMF) . 

In this paper, again using the same set of parameters, we show 
that drift can cause significant differences in the radial and latitu- 
dinal dependence of cosmic rays for consecutive solar minimum periods. 
We also searched the literature for additional modulation features 
related to the IMF polarity reversal and therefore relevant to deter- 
mining the role of particle drift in the modulation of cosmic rays. 

2. Model Calculations 


Using the parameters of Part 1, we calculated the radial and lati- 
tudinal dependence of the 100 MeV proton intensity for two consecutive 
solar minimum periods. The radial dependence, with the radial gradient 
G in %/AU, is shown m Fig. la. The latitudinal dependence, with the 
latitudinal gradient G m %/degree, is shown for 1 AU and 20 AU in 
Fig. lb and Fig. 1c Respectively. For the D(+) period (1976), G 

remains almost constant for most of the heliosphere with a rapia 
increase near the boundary, and G almost identical at 1 AU and 20 AU. 
For the D (-) period (1965), G i§ overall larger and decreases toward 
the boundary, while G differs r signif icantly for the two drift cases at 
1 AU, and at 20 AU. 0 (Note the difference in the intensity for the two 
drift solutions at 1 AU and 20 AU) . The general behavior of the 
density gradients displayed in Fig. la-c is a feature of all drift 
models. 

One of the best manifestations of the effect of drift becomes 
apparent when the transport equation is solved using ©-(polar angle) 
dependent Gaussian spectra at the outer boundary instead of the full ©- 
independent interstellar spectrum. By shifting these spectra (half- 
width of 10°) consecutively in 10° intervals from the heliospheric 
poles to the equator, the results shown in Fig. 2 were obtained. 
Protons reaching Earth from a particular region on the outer boundary 
are expressed as percentage of those that would have reached Earth if a 
full interstellar spectrum was used. The calculations were repeated 
for increasing values of the parallel diffusion coefficient K fl with 
other parameters unchanged. Solution (b) in Fig. 2 corresponds to the 
solutions m Fig. la-c. (Note that (K n ) is a constant in the expres- 
sion for K„ given by Potgieter and Moraal, 1985). These results show, 
in contrast with the no-dnft case, that drift causes positively 
charged particles reaching Earth to predominantly come from the outer 
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equatorial regions during a D(-) period (1965), but from the outer 
polar regions when the IMP changes polarity. A factor 3 increase in K„ 
extends this region to include mid-latitudes, while the situation 
remains virtually unchanged for the D(-) period. In general these 
calculations show that oppositely charged particles should traverse 
different regions of the heliosphere (Part 1). Another implication is 
that the (anti) correlation between variations in the CR intensity and 
solar activity parameters should exhibit a 22-year cycle. 



Polar Angle 9 (Degrees) Polar Angle 0 (Degrees) 


Fig. is (a) The radial dependence 
m the equatorial plane (0 = 90°), 
(b) the latitudinal dependence at 1 
AU and (c) at 20 AU, of the 100 MeV 
proton intensity. (The differen- 
tial gradients are m %/AU; -lati- 
tudinal gradients are in %/degree) . 


Fig. 2: The calculated percentage 
of 1 GeV protons reaching Earth 
from a particular region at the 
outer boundary for the two drift 
cases compared to the no-drift 
situation. The equatorial plane 
is at 0 = 90°. 


3. Discussion of Results and Observations 

Fig. la-c show that drift can cause significant differences in 
the radial and latitude dependence of cosmic rays during consecutive 
solar minimum periods. Compared to our results, the observed integral 
latitudinal gradient (Decker et al., 1984) is as yet inconclusive about 
the role of drift. Jokipn (1984) , however, by carefully simulating the 
exact observational conditions near the neutral sheet using a 
three-dimension drift model, found excellent agreement with the results 
of Newkirk and Fisk (1985) , who studied the statistical dependence of 
CR intensity on the distance from the neutral sheet. 

Less encouraging are observed integral radial gradients recently 
reported by Webber and Lockwood (1985) and McKibben et al. (1985). The 
radial gradient seems to decrease since 1981, which is not expected 
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from a drift point of view. It also appears that our model predicts 
the radial gradient too small for the 1976 solar minimum period. 
Since both the radial and latitudinal gradients are rather insensitive 
to parameter variations in drift models, these observations may become 
an interesting challenge to these models. It is also a complicating 
factor m explaining the role of drift concerning the anomalous compon- 
ents. However, the near constancy of the integral radial gradient 
observed before 1980, despite large changes in solar activity (Webber 
and Lockwood, 1985) , is consistent with our model, but not with 
time-independent spherically symmetric models. In these models a 
direct relation exists for variations in K n and the radial gradient. 

Potgieter et al. (1980) calculated neutron monitor differential 
response functions from latitude survey data for 1954, 1965 and 1976, 
and found that the 1965 response function deviates significantly from 
those for 1954 and 1976. Using the parameters of Part 1, which fit the 
low-energy 1965-66 and 1977 proton and electron data, we could not 
simulate the observed large difference between the corresponding 
response functions, although we obtained a split m the right direc- 
tion. The complicating factor here is that the low-energy proton 
fluxes were higher in 1976-77 than in 1965, while the neutron monitor 
counting rates were lower m 1976-77 compared to 1965. 

Another relevant observation is the change in direction of the 
annual wave vector, derived from the yearly variations in CR intensity, 
m 1958-59 and again in 1968-69 (e.g., Antonucci et al., 1978? Nosaka, 
et al., 1984). Whether drift effects are the predominant cause of this 
observation, is not yet clear and has to await more detailed studies. 

From 1972-77 the CR intensity was an extended, rather flat plateau 
preceded by a fast recovery over 2 years m 1970-72. (See the 
correlation study of Akasofu et al., 1985 regarding this period). 
During the 1965 solar minimum the intensity was peaked, preceded by an 
extended recovery period of ^ 7 years. This behavior is consistent 
with the feature of drift models that the proton intensity is less 
sensitive to variation m modulation parameters during a D(+) cycle 
(1970-80) compared to the 11-years before and after this period 
(Jokipn and Thomas, 1981? Kota and Jokipn, 1983; Potgieter and 
Moraal, 1985). 

Shea and Smart (1981) found a correlation coefficient of -0.86 
between the geomagnetic aa-mdex - a measure of disturbance m the 
ecliptic plane - and the CR intensity for the years around 1965. For 
this period Aldagarova et al. (1979) found that the CR intensity cor- 
relates best to coronal green line (CGL) intensity in the -20° to +20° 
heliolatitude range. However, for the years around 1976, Shea and 
Smart found a correlation coefficient of +0.28 and Aldagarova et al. 
the best correlation for +30° to +50°, and to a lesser extent for -20° 
to -30° heliolatitude. In combination these observations indicate that 
during a D (+) cycle (1976), variations in the CR intensity correlate 
better with solar activity parameters m a much wider heliolatitude 
range compared to a D(-) period (1965) when cosmic rays are predomin- 
antly transported to the inner heliosphere via the equatorial regions. 
(See also Jokipn, 1981). The same conclusion can be reached from the 
observations of Nagashima and Monshita (1980), Pandey et al. (1983) 
and Vernov et al. (1983) , despite different data and techniques used by 
them. These observations are consistent with our calculations shown m 
Fig. 2 and the drift velocities shown in Part 1. 
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4. Summary and Conclusion 

In Part 1 we discussed two key observations relevant to 
determining the relative importance of drift in cosmic-ray modulation. 
In the present paper, using the same set of parameters, we have 
illustrated the significant effect of drift on the radial and 
latitudinal dependence of cosmic rays for consecutive solar minimum 
periods. Compared with the integral radial gradient observed m 1976 
(Webber and Lockwood, 1985), the calculated value seems too small. A 
detailed comparison will however have to await the forthcoming solar 
minimum. The same applies to the latitudinal gradient which is as yet 
inconclusive about drift effects. 

Searching the literature for observations related to the IMF 
polarity reversal (in addition to Part 1) , we found distinct difference 
m neutron monitor response functions for consecutive solar minimum 
periods, and also m the annual variations of cosmic rays observed 
before and after polarity reversals. Whether drift is the predominant 
effect is however not yet clear. We also found several reports which 
indicate better correlation between variations m the cosmic-ray 
intensity and solar activity parameters (e.g., the corona green line 
intensity) over a much wider range of heliolatitude during 1970-80 
compared to before this period. These observations are consistent with 
drift models according to which cosmic ray protons primarily come via 
the polar regions during a D(+) period (1970-80), but primarily via the 
equatorial regions during D(-) periods. The observed peak vs. plateau 
in cosmic-ray intensity for the years around 1965 and 1976 respectively 
is also consistent with the general feature of drift models according 
to which the proton intensity is more sensitive to changes in 
modulation conditions during 1970-80 compared to before 1970. 
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